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I.  Foreward 


In  1962  we  made  a  preliminary  analysis  [2]  on  the  general¬ 
ized  aerodynamic  and  thermal  properties  of  hydrogen-fueled  high 
Mach  number  ramjet  engines  based  on  the  thermodynamic  properties 
and  computation  method  given  in  Reference  [1].  Since  more  de¬ 
tailed  characterization  [3]  became  available  on  the  thermal  pro¬ 
perties  of  hydrogen-air  combustion  products,  we  have  performed 
computation  checks  and  extended  our  calculation  to  the  effects 
of  non-equilibrium  flow  [4].  In  the  meantime,  a  master  curve 
method  was  developed  in  the  course  of  improving  our  calculation 
scheme  and  the  characteristics  of  equilibrium  and  nonequilibrium 
expansion  of  hydrocarbon  fuel  has  been  calculated.  (See  Part 
la  in  appendix). 

This  report  gives  a  composite  account  of  the  calculation 
results  mentioned  above  and  discusses  the  following  major  prob¬ 
lems  : 


(1)  Investigate  the  effects  on  the  aerodynamic  and  thermo¬ 
dynamic  characteristics  of  the  engine  by  various  working  para¬ 
meters  and  design  parameters  such  as  the  Mach  number,  the  recov¬ 
ery  coefficient  of  inlet  pressure,  the  type  of  fuel  -  hydrogen 

or  hydrocarbon,  the  degree  of  expansion,  the  nonequilibrium  expan¬ 
sion  of  high  temperature  combustion  gas,  and  so  on. 

(2)  The  effects  of  pressure  Increase  under  the  wing. 

(3)  Investigate  the  operating  parameters,  aerodynamic  and 
thermal  properties,  and  adjustment  range  of  various  major  comp¬ 
onents  of  the  engine  under  the  requirements  of  maximum  thermal 
efficiency  and  thrust  coefficient  while  keeping  the  design  rel¬ 
atively  simple  (e.g.  unadjustable  nozzle  throat). 

Based  on  the  analysis  and  calculation  of  the  above  aspects, 
this  paper  presents  an  evaluation  of  the  conclusions  found  in 
the  literature  dealing  with  the  performance  and  application  of 
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ramjet  engines.  We  also  present  some  preliminary  information 
on  key  questions  in  research  of  high  Mach  number  ramjet  engines 
and  its  aerodynamic  and  thermodynamic  properties. 

A  Cross-sectional  area  of  passage 
a  Sound  velocity,  meters/sec 
CF  Thrust  coefficient  (=  Ft/^W1Ae/2gc 
Cw  Velocity  coefficient 
f  Fuel  to  air  weight  ratio 
Ffc  Thrust,  Kg 
gc  9.81  Kgm/Nsec2 
H  Flight  altitude 

Same  as  J(-iRp),  heat  release  of  unit  mass  of  fuel 
I  Specific  impulse  of  air,  N/ (.Kg/sec) 

CL 

If  Specific  impulse  of  fuel,  N/ (. Kg/sec) 
i  enthalpy  Kcal/Kg 

J  Mechanical  equivalent  of  heat,  426.9  Nm/Kcal 
M  Mach  number 

2 

P  Pressure,  N/cm 
q  ^/W2/2gc 

s  Entropy  Kcal/Kg°K 
T  Absolute  temperature,  °K 
W  Relative  flow  velocity,  m/sec 
j3  Fuel  coefficient  or  fuel  equivalent  ratio 
5  Deflection  angle  of  air  flow  around  the  shock  wave  at 

the  front  of  the  wing,  with  engine  mounted  under  the  wing. 
X  Specific  heat  ratio 
^  back  pressure  recovery  coefficient 
7|c  combustion  efficiency 

Kinetic  energy  efficiency  of  inlet 
7|c  Total  efficiency  of  engine 

Subscripts  - ,  various  cross  sections  of  the  engine. 

See  (a)-(f)  in  Figure  27  of  Ref.  (1). 
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0  static  parameters 

*  critical  parameters 

c  combustion  chamber  of  complete  expansion 
d  inlet 

e  equilibrium 

f  expansion  under  fixed  constituents  and  equilibrium  vibra¬ 

tion 

ff  expansion  under  fixed  constituents  and  vibration 

i  incomplete  expansion 

s  isentropic  process 

II.  Original  Data  Used  in  Characteristics  Calculations 

( 1)  .  Flight  path 

The  great  majority  of  calculations  (throughout  Ref. [2]  and 
[4]  and  also  in  part  of  Ref.  [53)  assumed  a  flight  path  where  H 
is  equal  to  24,  31  and  26  Km  for  M  of  3,  5  and  7  respectively. 
Judging  from  the  data  compared  in  Figure  5  of  Ref  [1],  the  flight 
starts  from  some  H-M  point  already  reached  and  approaches  q^  =  0 
(illegible)  and  the  cruise  flight  path  suggested  by  [11]  in  Ref¬ 
erence  [1].  Part  of  the  calculations  in  Reference  [53  assumes 
a  flight  path  of  p^O  between  1  and  (illegible)  atmospheres  (see 
Figure  2  in  la). 

(2)  .  Inlet  state  1  and  1°9 


Table  1  gives  the  atmospheric  parameters  and  the  static 
state  for  Mach  numbers  3,  5  and  7  on  the  first  type  of  flights 
described  above.  The  atmospheric  parameters  were  based  on  old 
data  which  are  somewhat  different  from  these  in  I.  The  static 
state  is  based  on  curves  in  Reference  [6]  and  is  less  accurate 
compared  to  that  obtained  from  the  thermodynamic  table  in  I. 
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Table  1.  Inlet  State 
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(3)  Recovery  coefficient  of  Inlet  back  pressure 


The  o-  M.^  curves  found  In  the  literature  show  considerable 
differences  (see  Figure  11  In  I).  Based  on  the  experimental  data 
reported  in  Reference  [7]  and  [8],  we  compiled  the  curve  In  the 
figure(  oactual)  and  take  It  to  be  the  maximum  value  of  that 
can  be  obtained  to  date.  In  addition,  we  have  also  chosen  the 
upper  and  lower  limits  of  the  upper  limit  curve  is  from  the 
lsentropic,  variable  Y expansion  pressure  inlet  of  Reference  [9] 
and  the  lower  limit  curve  corresponds  to  the  double  cone  inlet  of 
Reference  [10].  (However,  for  <  4,  the  value  given  in 
Reference  [10]  seems  too  high,  we  therefore  use  the  =  0 

multiple  wave  inlet  data  of  Reference  [11]). 


In  order  to  reduce  the  amount  of  computation,  we  have  only 
computed(  a  actual)  for  the  case  of  pressure  increase  under  the 
wing.  For  this  case  (  oactual)  is  obtained  from  the  curves  des¬ 
cribed  above  and  the  Mach  number  behind  the  lower  shock  wave. 
Using  the  results  of  Reference  [12]  and  choosing  the  shock  wave 

4 


deflection  angle  of  5  =  10,  various  parameters  for  the  flow  pass¬ 
ing  the  shock  wave  at  different  M1  values  can  all  be  obtained  from 
the  curves  of  Reference  [6].  (The  effects  of  real  gas  properties 
of  air  such  as  variable  specific  heat  have  been  taken  into  ac¬ 
count  ) . 

(4)  Back  pressure  recovery  coefficient  of  the  combustion 

chamber  a 

-  c 

a  is  estimated  from  the  Mach  number  M~  =  0  •  15  at  the  com- 
c  2 

bustion  chamber  inlet  for  a  constant  diameter  chamber  with  temp¬ 
erature  increase  and  friction  loss  taken  into  consideration.  The 

values  of  o  are  listed  in  Table  2. 
c 

Table  2. 


0  -  6 

1  0 

2-0 

3 

5 
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7 

1  (  ! 
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1 


(5) 


Combustion  Chamber  Outlet  Back  Pressure 


Based  on  the  choices  of  data  and  the  procedure  just  des¬ 
cribed,  the  P^*/pto(  =  )  values  obtained  are  given  in  Tables 
3A  and  3B. 
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Table  3a.  a  values  without  pressure  increase  under  the  wing 
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Table  3b.  a  values  with  pressure  increase  under  the  wing 


M1 
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m 

0.5  1.0  2.0 
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(6)  Fuel 

Most  of  the  calculations  were  made  for  Hg  fuel  because  of 
the  following  reasons  (1)  this  fuel  provides  a  greater  specific 
thrust,  (2)  it  is  favorable  for  regenerative  cooling  under  high 
conditions,  (3)  the  flow  of  the  hot  combustion  products  in 
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the  tail  ejector  is  close  to  being  an  equilibrium  expansion 
process.  (IV)  Reference  [11]  has  made  rather  complete  analysis 
for  kerosene  fuel  where  as  similar  information  for  H2  fuel  has 
been  lacking,  (V),  In  Reference  [13]  -  [16],  calculations  made  on 
H2  fuel  ramjet  engines  of  different  characteristics  have  led  to 
quite  different  results  regarding  launch  payload  capacity  and 
economy.  In  order  to  evaluate  these  differing  conclusions,  one 
must  judge  upon  the  feasibility  of  the  various  characteristics 
of  the  H2  fuel  ramjet  engine  assumed  in  these  references.  In 
this  calculation,  the  heat  release  of  H2  is  taken  to  be  = 
J(-iRp)  =  12.1  x  10^;  and  in  the  calculation  for  (CH2)L  fuel, 

Hp  is  taken  to  be  4.419  x  10^  N-m/Kg. 


(7 )  Calculations  of  thrust  and  thermodynamics  properties 


The  following  equation  is  still  widely  used  In  current  refer¬ 
ences  for  the  calculation  for  supersonic  ramjet  engines. 


(l-QG^Vg 

Is 


So 


a  . 


(i) 


It  has  been  pointed  out  in  section  II. 5  of  I  that  the  above 
equation  represents  net  thrust  only  when  the  flow  outside  the 
airplane  body  is  isentropic.  This  condition  is  approximately 
satisfied  only  for  subsonic  and,  for  supersonic  flight,  there 
is  Invariably  a  shock  wave  system  outside  the  plane  body.  There¬ 
fore,  the  effective  thrust  under  supersonic  conditions  generally 
can  no  longer  be  found  from  Equation  (1)  but  Instead  is  dependent 
upon  the  actual  shape  of  the  airplane  body,  for  the  sake  of 
comparison  with  data  reported  In  the  literature  and  avoiding  the 
difficulties  associated  with  resistance  computation  for  a  part¬ 
icular  body  design,  we  will  consistently  resort  to  the  follow¬ 
ing  calculations: 


Without  using  pressure  increase  under  the  wing 
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_(l+f)G1\Ve  G1W1 


gr 


g„  ^ p«Ae_pi  A  t  ^ 


With  pressure  increase  under  the  wing 

(1+T)G  Y7  G  w  (3) 

?t,i - g - +  pPAq  -(—  +PjA,)  cos  s 

C  °c 

w 

(Derivatives  of  Eqs .  (2)  and  (3)  can  be  found  in  Section 
II. U  of  I). 


In  Eq.  (3),  parameters  for  cross  section  1  can  be  computed 
from  shock  waves  under  the  wing  and  it  is  also  assumed  that  the 
direction  of  the  exhaust  and  the  resulting  internal  thrust  are 
both  parallel  to  the  direction  of  the  flight. 


Equations  for  other  relevant  thermodynamic  properties  are: 

,  w  ,  ^ 

ia,i  -  M 


z^crAt/y\ 


no,i 


xr wi 


.  r  W  1  A„ 

cv  ±  i =  2[(1+f^"!r  ~  rr 
*  t1  -7— p  v;  a  1  1  '  e 


2gc!  lie 


'c  is  taken  to  be  1  in  all  calculations  here.  When  compar¬ 
isons  are  to  be  made  with  other  authors,  the  results  contained 
in  this  paper  need  to  be  converted  assuming  =  0.95  ana  cal¬ 
culating  the  thrust  according  to  Eq .  (1).  The  total  efficiency, 
as  defined  in  Reference  [11],  is  given  by 
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’ll  ■ 


rf  "'l 

+it 


therefore,  in  comparison  with  the  results  of  Reference  [11], 
we  convert  the  values  in  [11]  as  follows: 


V  -cn^ 


c  ± 


where  2  g^H^,  =  2x32.2x778x18630  =  934  x  10°. 

C 

We  obtained  the  conversion  coefficient  'J  for  flight  paths 

with  dynamic  pressure  /^'"f  lb  as  given  in  Reference  [11]. 

a  - - =350 — a 

2gc  ft* 

The  values  of  C^_  are  tabulated  below. 

Table  4.  Cvj  values  for  =  350  lb/ft? 
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1 

5 
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m 

3 

9 

1  o 
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1  .028  ; 

d.041 

i  .057 

1  .076 

1  .090  ] 

1.126 

( 8 )  Complete  expansion  and  Inc-  onplete  expansion 


We  use  the  following  definitions  for  complete  expansion  and 
incomplete  expansion. 


Complete  expansion  -  when  the  gas  pressure  Pe  over  the  cross 
section  of  the  tail  ejector  pipe  outlet  is  equal  to  the  ambient 
atmosphere  pressure  p^.  Actually,  in  supersonic  flights,  the 
back  pressure  behind  the  tail  ejector  pipe  is  not  equal  to  p-^ 
and  this  is  especially  true  when  there  is  a  pressure  increase 
under  the  wing.  For  this  reason,  the  definition  given  here  for 
complete  expansion  should  be  regarded  as  a  comparative  standard. 


Incomplete  expansion  -  consider  only  one  type  of  tail  eject¬ 
or,  the  ratio  of  the  outlet  cross  sectional  area  and  the  gas  col¬ 
lecting  area  Ae/A^  =  1.5*  The  ratio  is  appropriate  for  the  ex- 
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terior  design  and  resistance  where  the  under  wing  pressure 
increase  is  not  used.  It  is  also  convenient  for  comparing 
the  results  here  and  those  of  Reference  [11]. 

We  calculated  the  properties  for  complete  expansion  be¬ 
cause  generally  one  is  likely  to  achieve  complete  expansion 
using  the  body  and  wing  area  facing  the  wind  while  not  nec¬ 
essarily  increasing  the  external  resistance. 

( 9 )  Equilibrium  and  nonequilibrium  expansion  processes 

We  have  made  calculations  for  hydrogen  fuel  and  hydro¬ 
carbon  fuel  under  complete  expansion  for  three  expansion  pro¬ 
cesses:  expansion  process  where  both  the  components  and  vi¬ 

bration  are  kept  under  equilibrium,  expansion  where  components 
are  frozen  (from  3°°  and  up)  and  vibrations  are  in  equilibrium, 
and  finally  expansion  process  where  both  are  frozen.  (Methods 
of  computation  are  given  in  Section  II. 4  in  I). 

Over  the  upstream  of  the  ejector  throat,  flow  is  assumed 
to  be  isentropic.  Cw  =  0.97  in  calculating  the  outlet  nozzle 
velocity  We. 

(10)  Calculation  Results 

We  have  made  computations  for  20  design  points  for  the 
configuration  without  under  wing  pressurization.  Thermody¬ 
namic  properties  of  the  engine  and  major  aerodynamic  and  thermal 
parameters  for  various  cross  sectional  areas  are  summarized  in 
the  attached  tables . (Table  1  and  Table  2).  We  have  plotted 
these  results  against  /Q  to  observe  the  values  when  ^Jo.i  and 
Cp/L  reach  their ' respective  maximum  values  for  different  M^, 
see  Figures  1,  3  and  4.  The  variation  from  maximum  efficiency 
to  maximum  thrust  coefficient  is  obtained  for  different 
values  by  plotting  the  aerodynamic  thermal  parameter  and  thermo¬ 
dynamic  properties  just  obtained  versus  M^,  see  Figures  2  and  5. 
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III.  Effects  of  M1  and  0  on  aerodynamic  thermal  parameter 
and  thermodynamic  properties 

To  observe  the  and  effects  on  the  aerodynamic  thermal 
parameter  and  on  the  thermodynamic  characteristics  of  the  engine, 
we  first  plot  the  data  [5]  (in  attached  Table  1)  for  the  case 
of  no  pressurization  under  the  wing  versus  0  .  See  Figure  la- 
le.  It  can  be  seen  from  the  figure  that  for  small  0,  the  gas 
temperature  TO  increases  rapidly  with  increasing  0  .  (The  in- 

o 

crease  is  less  rapid  for  large  0).  The  rise  of  T^°  slows  down 
as  0  approaches  1  because  of  the  boundary  at  high  temperatures. 
T^O  reaches  its  maximum  in  the  vicinity  of  0  =  1.1,  while  T^O 
follows  the  increasing  0  all  along.  Under  these  conditions,  the 
increase  in  We  is  rather  small  (see  Figure  lb)  after  0  >1.28 
even  if  combustion  gas  can  achieve  complete  equilibrium  expan¬ 
sion  in  the  exhaust  nozzle.  In  certain  cases,  We  can  even  de¬ 
crease. 


For  the  same  0,  T^°  and  We  increase  with  M^,  since  T-^0 
increases  with  M^.  In  the  meantime,  both  (We  -  W-^)  and  VIe/W1 
are  decreasing  with  M-^  as  shown  in  Figure  lc .  As  a  result, 

I  ,  and  I„  .  decreases  with  M. ,  as  shown  in  Figures  Id  and  le. 

a-,1  i  -,1  l 

Figure  le  also  shows  that  If  ^  reaches  its  maximum  for  larger 
0  as  increases,  however  t?se  maximum  0  value  (  "0.6)  is  still 
considerably  less  than  1  even  for  =  7 . 


v, 

Since  (u,i  is  proportional  to  the  product  of  Ir  .  and  W  , 

Yl  1  1  1 

iiu  increases  with  increasing  and  for  a  given  M^,  the  0 
value  for  maximum  is  the  same  as  for  1^.  Calculations 

for  hydrogen  fuel  under  equilibrium  flow  and  Complete  expansion 
conditions  lead  to  the  following  approximate  0  values  at  maximum 
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,  depicted  in  Figure 
=  )  and  Af=£  2 
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The  value  of  Ae/A 

by  W,/W  (Figure  lc)  and 

x  e  re  *-  r<>  r  i  J 

effect  on  Ae/A-^  and  the  variation  with  is  similar  to  that  of 

We;  as  a  result,  the  changes  in  Cp  i  are  similar  to  those  of 


has  only  a  small 
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For  M-^  =  5  and  M1  =  7 


the 


Variations 


are  extremely  small 


for  g  greater  than  1.5  . 


Figure  li  shows  the  variation  of  yet  another  thrust  coef- 
ficient  F  (P]_A-^ )  whose  value  is  determined  mainly  by  and  We/ 
W^.  Therefore,  the  variation  of  Ft/P]_A2.  with  8  is  similar  to 
that  of  We/W^  for  the  same  and  increase  with  for  the  same 

0  . 


The  variation  [2]  of  A^/A2  as  a  function  of  is  shown  in 
Figure  lj .  Because  the  rate  of  increase  of  the  density  ratio 
(FiiP,)  is  much  faster  than  that  of  velocity  ratio  (W1/W2),  k^/k^ 
increases  rapidly  as  gets  larger.  It  seems  that,  for  an 
engine  with  constant  A2,  the  values  of  A^,  needs  to  be  increased 
7  times  as  M1  goes  from  3  to  7.  It  can  also  be  seen  from  Figures 
lk  and  11  that  yd  does  not  large  effects  on  A^#/A2  and  Ae/A2, 
however,  Ae  again  needs  to  be  increased  by  a  factor  of  7  as  M 
increases  from  3  to  7. 

IV .  Aerodynamic  and  thermodynamic  properties  at  maximum  ^°- 
and  maximum  for  different  M  numbers 


From  the  three  curves  for  equals  3,  5  and  7  in  Figure 
1,  we  have  chosen  Q.l,  0.3,  0.6  and  2.5  as  the  y3  values  for 
maximum  ^lo,^  and  maximum  Cp ,1%  and  indicated  these  in 

Figure  2.  Figure  2A  indicates  that,  as  increases,  the 

maximization  of  both  q ,  and  Cp.L,e  requires  T^O  to 

increase,  although  the  former  has  a  much  smaller  effect.  When 
M-^  is  slightly  greater  than  7,  T^O  for  both  are  the  same  even 
though  the  corresponding  fi  values  are  far  apart  (~*'0.6  and~2). 

This  is  because,  at  high  temperatures,  T^o  decreases  as  $  In¬ 
creases  (see  Figure  1A).  If  the  flow  can  be  kept  in  equili¬ 
brium,  the  i^O  value  at  high  T^O  can  still  increase  the  We 
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value.  Therefore,  when  M^  =  7,  the  difference  in  We  for  maxi¬ 
mum  *1  c.ifc  and  for  maximum  is  greater  than  that  of 

T^O.  From  Figure  2C,  we  can  see  that,  as  increases  from  3 
to  7,  the  value  of  We/Wp  for  maximum  is  almost  a 

constant  1.2  whereas  the  We/W^  for  maximum  has  de¬ 

creased  from  2.5  to  1.4.  This  variation  of  We/W..  basically 

t  1 

determines  the  variations  of  maximum  f.i.e  and  maximum 
7ja  with  Mp.  From  more  than  6000,  the  maximum  t.e  de” 

creases  by  a  factor  of  2  (Figure  2f)  while  maximum  ^lu-i  con¬ 
tinues  to  climb  with  and  reaches  its  peak  value  of  0.626 
near  M,  *  7  (Figure  2g). 

Since  the  Ae/Ap  for  maximum  Cpt-e  increases  very 
slightly  with  My ,  the  maximum  Cf./e.  ,  like  We/Wy,  decreases 
rapidly  with  M-,  (from  0.88  to  0.24).  For  high  M-^  values,  one 
effective  way  to  increase  Cp  is  to  use  incomplete  expansion  and 
thus  a  reduced  Ae . 

V.  Inlet  Stopping-Pressure  Recovery  Coefficient  and  Complete¬ 
ness  of  Expansion  -  their  Effects  on  the  Aerodynamic  and  Thermo¬ 
dynamic  Properties 

The  influence  of  high  ad  value  and  low  value  on  the 
thermodynamic  properties  can  be  realized  in  Figure  3  (computation 
of  Reference  [2])  and  Figure  6.  For  Mp  greater  than  5,  the  in¬ 
creasing  has  a  gradually  diminishing ' effect  on  tfo.i.e  and 
If  !_'  £  .  There  is  only  about  5%  increase  in  ^oA,e  or 

as  a  increases  from  a  actual  to  a  high.  (The  effect  is  greater 
under  nonequilibrium  conditions  at  high  temperatures.) 

Figures  3  and  5  show  the  effects  of  increasing  Cp,  decreasing 
nozzle  outlet  cross  sectional  area  Ae  on  the  completeness  of 
expansion.  Fixing  the  Ae  value  at  1.5A,  for  all  M^  values  leads 
to  the  following  results  (1)  increase  of  Ae  with  is  very  small 
(Figure  3d),  (2)  We  decreases,  the  amount  of  decrease  is  smaller 
for  high  than  for  low  (Figure  3A),  (3)  increases, 
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the  amount  of  increase  is  smaller  at  higher  M^,  see  Figure  3c, 
(4)  ^lo,c,e  decreases,  the  amount  of  decrease  is  greater 

for  higher  except  when  3  >  1.5*  We  noticed  that,  when 

=  7  and  ji  =  0.5,  setting  A&  equal  to  1.5A,  will  decrease 
*1*1*  by  282 * 

VI.  Effects  of  Fuel  and  Nonequilibrium  Expansion 


Although  the  discussion  so  far  has  been  primarily  based  on 
the  calculation  results  for  hydrogen  fuel,Q^4j  the  data  for 
hydrocarbon  fuel  under  the  same  flight  conditions  B  have  also 
been  included  in  Figure  1.  The  following  observations  can  be 
made  from  the  graph. 

Cl)  For  the  same  T-^0  and  3,  the  T^O  values  for  the  two 

types  of  fuel  are  very  close  to  being  equal  even  when  the  fuel 

amounts  differ  by  more  than  a  factor  of  2  . 

(2)  For  the  same  T.^0  and  3,  (i^O  -  ig )  for  hydrogen  fuel 

is  greater  than  that  for  hydrocarbon  fuel  and  the  difference  gets 
larger  at  higher  temperatures.  Thus,  the  equilibrium  expansion 
ejection  speed  of  hydrogen  fuel  is  greater  compared  to  hydro¬ 
carbon  fuel,  and  the  difference  in  ejection  speed  Increases  as 

or  T^O  increases  (Figures  IB  and  1C). 

(3)  The  low  I  value  of  hydrocarbon  is  due  to  its  lower 
We,  see  Figure  Id;  however,  because  hydrocarbon  has  a  lower  fH^, 
than  hydrogen,  the  result  is  that  the  overall  efficiencies  g 
for  both  fuel  are  quite  close  (Figure  If). 

(.4)  Because  of  the  lower  I  and  higher  (more  than  twice) 
f  of  hydrocarbon  fuel,  its  is  lower  than  that  of  hydro¬ 

gen  by  more  than  a  factor  of  2  for  the  same  and  3  .  (Figure 
le) . 

(5)  For  the  same  T^O  and  P  ,  the  molecular  weight  of  hydro¬ 
carbon  fuel  gas  is  higher  than  the  molecular  weight  of  hydrogen 
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by  approximately  20%  (Figure  7).  Thus,  for  the  same  Tg  and  Pg 
the  density  of  the  former  is  20%  higher  than  the  latter.  The 
ejector  cross-sectional  areas  of  the  former  are  2055  smaller 
than  those  of  the  latter  for  the  same  flow  rate  and  ejection 
speed.  Judging  from  Figure  lg,  the  ejector  cross  section  for 
hydrocarbon  fuel  is  approximately  1055  smaller  than  its  counter¬ 
part  for  hydrogen  for  the  same  M1  and /3 . 

Perhaps  another  major  difference  between  the  two  types  of 
fuels  can  be  seen  from  the  effects  of  nonequilibrium  flow.  Fig, 
8B  Indicates  the  effects  of  nonequilibrium  expansion  on  the  char¬ 
acteristics  of  the  two  fuels  upon  approaching  maximum 
under  various  conditions.  For  high  and  high  /}  ,  the  per¬ 
centage  decrease  in  Wg  of  hydrocarbon  due  to  nonequilibrium 
flow  is  not  much  larger  than  por  hydrogen  fuel,  nevertheless, 
this  decrease  is  weighted  heavier  for  Wg  -  W, ,  and,  as  a  result, 
the  decreases  in  Iai>  If  f,  *lo.L  are  larger  for  hydrocarbons. 

See  Figure  8d-8f.  1 

VII.  Influence  of  altitude 


As  an  example  in  I  ,  we  obtained  the  aerodynamic  and  ther- 

CL 

modynamlc  properties  of  the  engine  when  is  kept  between  1 

and  10  atm..  See  Figures  8-10  In  la.  The  equilibrium  and  frozen- 
constituent  and  If  versus  under  maximum  Tly.  and 
maximum  ,  as  shown  in  Figure  9,  can  be  obtained  from  the 

data  in  Figures  8-10  using  the  envelope  curve  of  different 76  '  s 
or  plotting  %  and  against  ft  .  Figure  9  clearly  shows  the 
advantages  of  lower  flight  altitudes  at  high  M1  —  not  only 

»  and  *4,6  are  increased,  the  effects  of  nonequilibrium 
flow  are  also  substantially  reduced.  It  seems  that  reducing 
the  flight  altitude  is  an  effective  means  of  minimizing  the 
nonequilibrium  effects  at  high  M^. 


VIII.  Comparisons  of  thermodynamic  properties  of  partially 
adjustable  and  totally  adjustable  engines  and  their 
adjustment  criteria. 

Since  all  our  calculations  are  made  under  the  assumption 
that  components  achieve  their  optimum  performance  for  a  given 
M^,  these  calculations  should  thus  be  viewed  as  corresponding 
to  engine  performance  where  the  geometric  configuration  is  tot¬ 
ally  adjustable  to  meet  optimization  criteria  when  there  is 
a  variation  in  M,  .  In  this  section,  we  will  make  some  analy¬ 
sis  and  comparisons  of  the  component  adjustment  criteria  and 
their  corresponding  thermodynamic  performances  for  the  follow¬ 
ing  three  situations:  (a)  requiring  the  maximum  efficiency 
for  all  flight  Mach  numbers  M^,  (b)  requiring  the 
maximum  thrust  coefficient  for  all  M-,,  and  (c)  allow¬ 

ing  a  constant  cross-sectional  area  at  the  ejector  throat. 

(1)  Maintaining  maximum  efficiency 

First,  the  maximum  efficiency  /3  -  M]_  curve  (the  ^ \ ‘  curve 

in  Figure  5a)  is  obtained  from  the  curves  for  various 

M^  under  complete  expansion  in  Figure  3.  The  variation  pattern 
of  ,  Cf..\e  ,  A^VA^  T^O,  and  for  maximum 

efficiency  can  be  found  from  other  curves  in  Figure  3*  Results 
are  shown  as  thick  solid  curves  in  Figure  5. 

(2)  Maintaining  maximum  thrust  coefficient 

Theyj  -  M^  curve  corresponding  to  maximum  Cp.i.e  is  ob¬ 
tained  from  the  curves  for  different  in  Figure  3 

where  Incomplete  expansion  with  Ag/A^  ■  1.5  is  assumed.  The 
result  is  the  line  in  Figure  5c.  The  variations  of  other 

parameters  under  the  maximum  Cp  condition  can  then  be  obtained 
from  other  curves  in  Figure  3.  See  thick  dashed  line  in  Figure 
5. 

(3)  Keeping  A^*  constant 


16 


The  thick  dash-and-dot  curve  in  Figure  5  shows  the  variation 
of  parameters  by  using  Figure  3  and  setting  A^*  *  0.34  . 

In  view  of  the  variation  pattern  of  various  parameters  un¬ 
der  the  three  adjustment  conditions  listed  above,  the  constant 
ejector  throat  scheme  is  only  slightly  behind  the  maximum  ef¬ 
ficiency  scheme  in  its  effectiveness,  thrust  performance,  T^O 
and  Ag  while  maintaining  its  advantage  of  not  having  to  adjust 
the  throat.  This  comes  about  because  the  maximum  efficiency 
scheme  does  not  require  any  large  changes  in  A^#.  (see  Figure 

5f ) .  As  compared  to  scheme  (2),  scheme  (3)  has  high  efficiency 
(upper  70"),  is  easy  to  adjust  and  its  combustion  temperature  T^O 
is  also  low.  A  low  T^O  implies  a  relatively  small  boundary  effect 
in  real  flow  and  less  problems  with  heat  transfer  and  material  stru 
ture.  Disadvantages  are  a  lower  Cp  ^  and  relatively  large  adjust¬ 
ments  in  outlet  cros sjsection  e. 

IX.  Thermodynamic  performance  and  adjustment  criteria  of  engines 
with  under  wing  pressurization. 

Table  2  (attached)  and  Figure  4  are  compilations  of  calcula¬ 
tion  results  when  under  wing  pressurization  is  used.  Based  on 
these  data,  the  engine  performance  under  the  three  adjustment 
conditions  discussed  above  can  be  arrived  at  in  a  similar  fash¬ 
ion  to  the  case  without  under  wing  pressurization.  These  re¬ 
sults  are  included  in  Figure  5  as  five  lines  [Note  -  the  compu¬ 
tation  of  Cf  i  with  pressurization  is  based  on  the  complete 
expansion  assumption]. 

(1)  Effects  of  underwing  pressurization 


The  advantages,  as  can  be  seen  From  Figure  5,  are 

(a)  The  increase  in  is  larse-  The  amount  of  inc¬ 

rease  and  the  percentage  increase  are  both  increasing  with  M, 

± 

(Figure  5b).  For  example,  for  M  ■  7 ,  the  efficiency  gained 


through  pressurization  is  equal  to 

*|a, =  (  0.19  ••  0.4»i?)-r  0-6IT  *  o.2S  s. 

The  increase  in  efficiency  is  probably  caused  by  two  fac¬ 
tors.  First,  the  total  tf-value  is  increased  because  of  the  ad¬ 
ditional  shock  wave  from  the  pressurization.  Secondly,  the 
resistance  produced  by  the  first  wave  is  Included  in  that  of 
the  wing.  Since  there  are  always  shock  waves  produced  by  the 
wing,  the  net  effect  is  equivalent  to  circumventing  the  resis¬ 
tance  of  the  first  wave. 

(b)  Maximum  Cp,t,€  (complete  expansion)  has  a  slight 
Increase,  see  Figure  5c. 

(c)  There  is  little  variation  in  A^/A^  This  is  more 
evident  at  higher  (see  Figure  5e).  As  a  result,  the  adjust¬ 
ment  of  inlet  is  relatively  straightforward  at  high  M  numbers. 

(d)  If  one  uses  the  ^1.0,4  *  ^|o(ma>r  scheme,  the  variation 

in  exhaust  cross  section  Ag  is  very  small,  see  Figure  5e. 

On  the  other  hand,  pressurization  also  brings  about  some 
disadvantages : 

(a)  If  the  no,  scheme  is  used,  the  varia¬ 

tion  required  in  the  throat  cross  section  is  larger,  see  Figure 
5f. 

(b)  If  one  uses  the  A^*  *  constant  scheme,  no  problems 

are  encountered  for  M^<  6  where  adjustments  are  simple  and 
is  high.  For  >  6,  however,  the  decrease  of  is 

rapid.  The  A^*  *  constant  scheme  is  therefore  unsuitable  for 

>  6.  Without  pressurization,  this  scheme  is  acceptable  be¬ 
cause  of  its  satisfactory  at  high  M^. 


13 


An  overall  view  indicates  that  pressurization  under  the 
wing  does  not  change  the  performance  parameters  in  any  quali¬ 
tative  way,  its  effect  is  mainly  quantitative  and  usually  leads 
to  improvements.  Generally  speaking,  the  advantages  are  more 
prominent  at  higher  M^.  Naturally,  in  order  to  realize  the 
advantages  in  practice,  considerations  must  be  given  to  the 
design  configuration  of  the  engine,  wings  and  body.  For  in¬ 
stance,  since  the  angle  is  smaller  at  higher  M^,  the  engine 
inlet  must  be  placed  farther  back  from  the  front  edges  of  the 
wings.  Detailed  analysis  on  this  problem  can  be  found  in 
Reference  [12]. 

(2)  Comparison  of  the  three  adjustment  schemes 

In  the  =  constant  scheme,  the  adjustments  required  are 
small  and  thermodynamic  properties  are  good.  decreases 

rapidly  for  greater  than  6. 

In  the  maximum  scheme,  A^*  is  required  to  vary 

over  a  large  range  —  A^*  must  be  cut  in  half  when  increases 
from  5  to  7.  Compared  to  the  other  two  schemes,  T^O  is  much 
higher  and  this  is  true  especially  for  low  . 

In  the  maximum  efficiency  scheme,  the  required  changes  in 
Ag  is  small  and  T^O  is  also  low  for  less  than  7  (more  prom¬ 
inent  at  lower  M-^).  These  features  are  favorable  in  minimizing 
freezing  loss  and  in  structured  strength  considerations.  In 

this  scheme,  the  changes  in  A,*  is  much  less  than  that  required 

C  J 

in  the  maximum  method.  Another  advantage  is  that  Ag  is 

almost  constant.  Ag  can  therefore  be  unadjustable  because  small 
variations  in  expansion  have  very  little  effect  on  the  perfor¬ 
mance  when  total  expansion  Is  approached.  From  the  adjustment 
point  of  view,  this  scheme  is  easy  to  satisfy. 
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X.  Comparison  of  subsonic  combustion  and  supersonic  combustion 

The  performance  data  of  supersonic  combustion  engine  cited 
in  References  [10]  and  [11]  are  quite  close  together,  as  shown 
in  Figure  10.  When  one  compares  their  data  with  our  calculated 
results  on  subsonic  combustion  (thick  line  in  Figure  10),  one 
can  make  the  following  observations:  For  equilibrium  flow  and 
less  than  8,  a  highly  feasible  subsonic  combustion  ramjet 
engine  is  superior  to  the  supersonic  combustion  case.  This 
superior  tendency  is  estimated  to  hold  possibly  to  M-^  =  10  if 
one  extrapolates  the  curve  in  Figure  10.  Even  under  the  least 
favorable  subsonic  combustion  conditions,  as  the  lower  six 
curves  in  Reference  [19]  show,  supersonic  combustion  is  still 
inferior  for  M.^6. 

Naturally,  for  the  same  altitude  the  subsonic  combustion 
suffers  a  higher  loss  to  the  real  nonequilibrium  flow  than  the 
supersonic  combustion  does.  But,  as  one  can  see  from  the  figure, 
the  nonequilibrium  flow  loss  of  a  high  subsonic  combustion 
using  hydrogen  fuel  is  much  less  than  the  loss  due  to  use  of 
hydrocarbon  fuel.  The  total  efficiency  is  approximately  equal 
to  the  supersonic  one  at  equal  to  7.  At  lower  altitudes, 
of  subsonic  combustion  is  higher  and  the  nonequilibrium 
flow  loss  of  hydrocarbon  fuel  is  quite  small  at  =  7.  It  is 
believed  that  subsonic  combustion  ramjets  can  be  used  at  lower 
altitudes  for  M1  above  5. 


XI .  Feasibility  evaluation  and  comparison  analysis  of  results 
in  current  references 

We  compiled  the  relevant  data  on  efficiency  and  specific 
thrust  1^  and  plotted  them  in  Figs.  11  and  12.  In  these  figures, 
the  numbers  assigned  to  each  curve  are  the  same  as  the  reference 
numbers  in  I.  Our  computation  results  are  plotted  as  heavy  lines 
in  the  figures.  We  have  done  so  to  facilitate  the  comparison 
and  evaluation  of  the  ramjet  engine  properties  found  in  many 


current  references.  The  feasibility  of  using  the  ramjet  eng¬ 
ine  as  a  satellite  launcher  and  its  ability  can  also  be  analyzed 
based  on  these  data. 

Regarding  the  performance  of  hydrogen  fuel,  our  calculated 
equilibrium  flow  value  is  somewhat  higher  than  that  given  by 
Reference  £20]  in  I.  All  the  results  given  in  Reference  [13] 
in  I,  in  the  =  3  to  7  range,  are  higher  than  our  calculated 
value.  As  for  the  data  in  Reference  [10]  of  I,  since  the  same 
nonequilihrium  flow  loss  has  been  used  for  all  the  perfor¬ 

mance  is  too  low  at  low  the  data  for  =  7,  however,  seem 
reasonable.  The  performance  data  at  low  Mj  as  reported  in  Ref¬ 
erence  [11]  of  I  are  even  lower.  Data  in  Reference  [8]  of  I 
are  too  low  over  the  entire  range  and  the  variation  trend  in 
Reference  [12]  is  incorrect. 

If  If  value  of  hydrocarbon  fuel  reported  in  Reference  [19] 
of  I  for  high  is  very  close  to  our  calculation  and  those  re¬ 
ported  in  [10]  and  [11]  are  somewhat  too  low. 

XII.  Concluding  remarks 

Based  on  analysis  of  the  calculation  results  presented  in 
previous  sections,  we  can  make  the  following  preliminary  deduc¬ 
tions  : 

Cl)  When  the  M  number  Increases  from  3  to  7,  thermodynamic 
properties  of  the  ramjet  engine  are  expected  to  improve  further. 
However,  when  M  is  greater  than  5.5,  proper  design  should  be 
made  to  avoid  increasing  losses  due  to  nonequilibrium  expansion 
and  boundary  combustion  gas. 

C2)  When  the  M  number  reaches  7,  the  intake  gas  temperature 
of  the  combustion  chamber  is  2300°K  and  the  maximum  efficiency 
and  maximum  thrust  require  the  exhaust  temperature  to  be  3000°K. 
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At  this  temperature,  the  boundary  effect  of  the  combustion  gas 
increases  with  decreasing  pressure.  At  36  km  altitude,  c  is 
approximately  4.4  atm.  and  the  nonequilibrium  expansion  loss 
of  hydrocarbon  fuel  is  likely  to  make  the  efficiency  at  =  7 
less  than  the  efficiency  at  a  6.  If  the  altitude  is  reduced 
to  29  km,  then  ^0  =  10  atm  and  the  nonequilibrium  expansion 

loss  is  greatly  reduced  and  the  efficiency  continues  to  increase 
from  =  6  to  =  7  and  beyond.  This  effect  is  more  effective 
than  changing  to  hydrogen  fuel.  (The  flight  altitude  can  be 
increased  to  32  km  if  under  wing  pressurization  is  used,  see 
Figure  4  in  la). 

(3)  It  is  not  necessary  to  use  supersonic  combustion  be¬ 
fore  M,  reaches  9  or  10,  hence,  the  research  work  in  this  area 

j. 

can  be  delayed  until  a  few  years  later. 

(.4)  The  advantages  of  hydrogen  fuel  over  the  hydrocarbon 
fuel  are  the  following:  hydrogen  fuel  specific  thrust  is  twice 
as  high;  although  the  efficiencies  under  equilibrium  expansion 
are  not  much  different,  the  efficiency  of  hydrogen  fuel  is  much 
higher  under  actual  expansion  condition;  and  finally,  hydrogen 
fuel  is  a  better  coolant. 

(5)  underwing  pressurization  allows  the  inlet  channel  to 
work  at  a  lower  and  less-variable  M  number,  more  favorable  for 
better  performance. 

(6)  without  under  wing  pressurization,  a  fixed  ejector 
throat  cross  sectional  area  under  different  flight  speeds  has 
a  small  effect  on  the  performance.  With  pressurization,  the 
throat  area  must  be  varied  according  to  flying  speed  in  order 
to  maintain  the  maximum  efficiency  over  a  wide  range. 

(7)  The  super  high  speed  ramjet  engine  has  been  proposed 
as  the  second  stage  for  a  satellite  launch  vehicle.  Its  carry¬ 
ing  capacity  has  been  estimated  in  the  current  literature.  Judg¬ 
ing  from  the  thermodynamic  data  of  the  ramjet  engine,  the  esti- 
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maces  in  References  [5],  [6],  [10]  and  [13]  (in  Part  I)  are  all 
considered  to  be  high.  Those  in  [8]  and  [12]  are  too  low  and 
estimates  in  [9],  [10],  [19]  and  [20]  are  the  most  reasonable. 

(8)  Master  curves  in  la  can  be  used  to  greatly  reduce 
the  amount  of  aerodynamic  and  thermodynamic  analyses  similar 
to  those  presented  here,  and  in  the  meantime,  increase  the  accu¬ 
racy  of  the  analysis.  These  curves  are  even  more  useful  for 
those  fuels  for  which  thermodynamic  properties  are  not  yet  avail¬ 
able  . 
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Table  1.  Aerodynamic  and  Thermodynamic  properties  of  hydrogen 
fuel  under  complete  expansion  without  under  wing  pressurization. 
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Supplemental  Information  to  the  article  "Aerodynamic 
and  Thermodynamic  Characteristics  of  Supersonic  Ramjet  Engines 
with  Mach  Number  between  3  and  7". 

Part  la 


(1)  For  figures  13(a)  to  13  (d),  the  temperature  of  the  CnH2n 
fuel  is  T^O.  This  figure  is  based  on  the  thermodynamic  prop¬ 
erties  given  in  Reference  [1]  and  the  combustion  gas  tempera¬ 
ture  is  calculated  for  T-^0  =  200Q°K  and  M  =  5  to  6. 


(2)  Figure  5(c)  gives  the 


value  for  equilibrium 


isentropic  expansion  of  the  combustion  gas.  Although  the  veloci¬ 
ty  ratio  of  a  boundary  combustion  gas  undergoing  equilibrium 
expansion  cannot  be  found  from,  the  kg  formula  like  the  variable 
specific  heat  combustion  gas  without  a  boundary,  it  can  be 
seen  from  Figure  5(c)  nonetheless  that  temperature,  pressure 
and  0  have  rather  small  effects  on  the  velocity  ratio  under  the 
same  KQ  e  and  the  same  expansion  pressure  ratio  *  c r/p  .  The 
heavy  curve  in  the  figure  repr-  sents  the  velocity  ratio  for  a 
given  Kq  g  and  under  different  temperature,  pressure  and  . 

For  Kq  e  ^*1.16,  the  error  is  less  than  2%. 


(3)  Figure  14  shows  the  isentropic  index  K  Q  of  the  combus- 

o  > e 

tion  gas.  It  Is  computed  using  Eq.  (15)  in  I  and  the  Ae  value 
in  (.1)  (Substitute  for  Figure  23  in  I). 


(4)  Figure  15  shows  the  molecular  weight  It  of  the  combustion 

gas.  Data  were  directly  taken  from  Reference  [1]  and  can  be 

used  in  computing  An  and  W.  The  Rvalue  can  be  read  to  0.02 

u  9q 

from  the  graph  and  this  is  sufficiently  accurate  for  the  above 
computations . 


The  procedures  of  using  the  above  graphs  are  as  follows: 

(1)  For  given  M  and  H  (height)  values,  p^o  is  obtained 
from  Figure  1  and  T10  can  be  found  from  Figure  2. 
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(2)  P^O  is  found  from  p^O  and  Q  values  of  the  inlet 
and  combustion  chamber. 

(3)  T^O  can  be  obtained  fairly  accurately  from  Figure  13 
for  a  given  temperature  of  T  2000°K 


(4)  Ko,e  can  be  found  from  Figure  14  using  T^O,  1^°  , 
and  • 

is  found  from  Figure  3(c)  using  Ko,e  and  the 
specified  expansion  pressure  ratio  Pq/Pj  of  the  ejector. 
[(W/A0)-j  ^  can  be  found  in  existing  function  tables], 

(6)  JU.  is  obtained  from  Figure  15  using  P  ,  fj*  ,  and  T30, 
First,  compute  Ag2  using  the  following  equation 

gJ-K**  *9-8<  *847-8  «  7V 

and  then  calculate  W^. 


All  the  graphs  give  parameters  of  CnH2n  combustion  gas 
for  the  following  values:  ^  =  0.25,  0.5,  0.8,  1,  1.5,  and  2; 
f^o  =  CO. 3),  1,  5  and  10;  and  T^O  =  1600  -2800°K.  In  case 
the  (  or  ^o.)  value  in  the  computation  is  not  listed  in  the 
graph,  then  one  can  read  off  values  for  the  ^  (  or  ) 

value  in  the  graph  and  plot  them  against  or  P^O),  the  quan¬ 
tity  value  is  then  read  off  for  the  desired  or  p 
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Appendix:  Calculation  of  the  combustion  gas  state 
near  the  ejector  throat 

In  ejector  tube  calculations,  one  has  to  find  the  critical 
parameters  at  the  ejector  throat.  This  calculation  is  a  ted¬ 
ious  iteration  process  for  decomposing  fuel  gas.  We  now  des¬ 
cribe  a  very  accurate  method,  using  the  existing  thermal  prop¬ 
erties  tables,  in  finding  the  p*  and  T*  parameters  at  M  *  1  for 
given  total  pressure  pQ,  total  temperature  TQ  and  .  (Further 
calculations  of  throat  area  will  then  be  possibl^. 

First,  choose  three  p  values  from  the  table  with 
in  the  vicinity  of  0.5  and  carry  out  precise  isentropic  cal¬ 
culations  respectively  by  three-point  interpolation  using  S  * 

S  and  tabulated  T  values, 
o 

This  leads  to  calculated  T,  i,  w  and  a.  Using  m  and  the 
three  p  values,  find  p*  and  T*  for  M  *  1.  Obtain  i»,  w#,  and 

a*  from  the  table  and  check  to  see  if  W#/a#  is  equal  to  1.  Gen¬ 
erally,  one  cycle  of  such  calculation  leads  to  a  discrepancy  be¬ 
tween  M  and  1  of  less  than  0.001. 


Actual  computations  have  indicated  that  the  computation 
method  given  in  Reference  [21]  for  finding  P#  and  T*  for  decom¬ 
posing  fuel  gas  really  is  not  much  more  accurate  than  the  follow¬ 
ing  simple  method:  using  given  ^3  ,  TQ  and  ,  find  Ke,0  using 

the  thermal  properties  table  directly  or  using  Figure  14  of 
this  work,  then,  the  following  formulas  allow  one  to  compute 
P»  and  T#  Xe.o 


\*s 

/  J*.  /  arS  Senerally  within  2%. 


Further  calcula¬ 


tion  on  throat  area  A#  may  have  a  greater  error  (5%  or  more)  be¬ 
cause  of  the  accumulation  of  errors. 
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